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Abstract 



ly-v ■ We derive a general expression for the expectation value of the phase acquired by a time de- 

pendent wave function in a multi component system, as excursions are made in its coordinate 
K^ \ space. We then obtain the mean phase for the (linear dynamic E ® e) Jahn- Teller situation in an 

in 

0^ ' electronically degenerate system. We interpret the phase-change as an observable measure of the 



effective nodal structure of the wave function. 
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In a recent publication a geometric (or Berry) phase was calculated for the wave function 
of a mult i- component closed system [l| . This differs from the usually considered situations, 
in which the Berry phases emerge from the wave-function of the system during the cyclic 
evolution of some external parameter. It is of interest to point out that in a famous prototype 
of a multi-component, closed system (an electronically doubly degenerate molecule) dynamic 
solutions for the E ® e linear Jahn- Teller effect (DJTE) were fully obtained as long ago as 
1957 j2, y]. The two parts of system were the electronic and ionic constituents of the 
molecule. Though this has received, as just noted, a complete treatment early on, the 
dynamic problem has not left the scientific agenda ever since. Descriptions of some of the 
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early refinements are found in two books |^, |5[ ; the most recent publication known to us 
and involving a variational treatment of the problem is in '^. The physical consequences of 
the Berry phase on the DJTE were clearly brought out by the late Frank Ham [7| and more 
recently in J8| , both of which papers showed (albeit under different physical conditions) that 
the value of the Berry phase may be critical in determining the order of energy levels in 
the closed molecular system. This phase is thus clearlv observable by experiment \j\. Its 
physical interpretation, essentially along the lines of \Bn la] , will be given later in this work. 

In |jy] an operator was proposed for the phase change (called "quantized phase") in a 
closed system. Here we shall derive an expression for this phase from first principles and use 
it to calculate the phase change in the vibronic doublet ground states oi an E ® e linearly 
coupled Jahn- Teller system. We shall use a "guessed solution", for the ground state, which 
is transparent, intuitively simple and algebraically easily manageable y, y, lll|. Though 
not variationally obtained, the "guessed solution" was found to have eigen-energies that 
are considerably closer to the exact, computed energies of [3| than any other approximate 
solution with which it was compared. This comparison is seen in Fig. 2 of [11]. Later 
treatments did not test their methods by comparison with the "guessed solution". 

Our point of departure is the time-dependent Schrodinger equation 

i-V'(r,t) = ff(r,t)^(r,t) (1) 

for a wave function '?/'(r, t) that depends on the internal coordinates r of the system, as well 
as on time t. {h = 1) The system is coupled to the environment; hence the dependence on 
t in the Hamiltonian. In a closed system, the Hamiltonian is time-independent, ip is still 
time-dependent, as, e.g., in a wave-packet. 



We write the state (assumed to be regular in the coordinate space and vanishing at its 
boundaries) as 

^ = Ae'^ (2) 

with A and S being real functions of r, t. {A > 0.) We shall utilize the equation of continuity 
and the Hamilton- Jacobi equation |15i|: 

±^^Sr^-?^-Vir,t)^±-A-WA (4) 

Here ?7i is a mass parameter common to all degrees of freedom, with all coordinates scaled 
to this mass. V is the potential. 

Let us now consider the change in the wave function, between the initial state of the 
system at t = and a final time tf. Real and imaginary parts of the change in the logarithm 
are 

[M*^ = [InMl' + t[arg{^)]l' = [lnA]l' + t[S]l' (5) 

These are functions of the coordinates. To form quantum mechanical expectation values 
(denoted by angular brackets about the relevant quantities) we multiply by A"^ and integrate 
over all coordinates. Thus the mean of the changes can be written as 

< [lnA\l^ >+i< [S]l'' >= I ^ dt— I dvA^ilnA + iS) (6) 

[ It is natural to conjecture that to form statistical expectation values (appropriate to mixed 
states |l2]), the factor A^ is to be multiplied by the relative statistical weight of the state 
and the contribution due to all states be summed over. We do not pursue this topic here.] 
Separating real and imaginary parts we get for the real part the quantity — |[5e]o''^ where 

S, = - j dr\^\MM'' (7) 

Se (which is different from the phase S) is reminiscent of a von Neumann entropy, in which 
the density operator is projected onto the initial state. 

We turn now to the rate of change of the expectation value of the phase ^^^- We change 
the order of integrations in equation (jH)) and obtain after some manipulations : 

15_^ = I dv{A'^ + S^)= fdr[A'^ - -SV{A'VS)] 
dt J ^ dt dt ' J ^ dt m ^ '^ 

dAA'^-^^-A\^Sf\ (8) 

at m 



having used equation Q and integrated by parts (with vanishing integrands at space- 
extremities). We now substitute for (VS*)^ from equation (0} and obtain a change in the 
sign of the first term in the above expression, as well as the expectation values of (twice) 
the potential and a term related to the kinetic energies, which can be reworked by a further 
integration by parts so as to put it into a form of definite sign, giving 

We next recall that S = Im{lnip) and reinstate the time integration to get the change of 
phase as 
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[<S>]y = - dt /dr[/m(V.+(r,t)^^) + 2|(^(r,t)|V(r,t) 



+ l(V|^(r,t)|)2] (10) 

m 

where the cross means Hermitian conjugate. On the other hand, multiplying equation (^ 
by ?/'"'" (r,t), integrating over the coordinates and again integrating by parts, we obtain 

/m|rfr(V^+(r,t)g^) = -|c^r '^^^^^^^'' -|cir|(V^(r,t)|V(r,t) (11) 

which we use to eliminate the expectation value of the potential from equation (|TO|l . We 
then get: 



[< S >\^ = ImJ dtj rfr(^+(r,t)^^^) 

+ 1 f' dt /rfr[|V^(r,t)r-(V|^(r,t)|)2] (12) 

m Jo J 

This is our central result for the mean phase change. The first term is of the form familiar 
in e.g., expressions of the open path geometric phase llSf. In the second term, to be denoted 
for brevity 6K, one has the difference between two space-derivative terms, one involving the 
total (complex) time-dependent wave-function and the other its modulus. 

We now calculate the phase change in our molecular model for a multi-component closed 
system. The main simplification in the model is the restriction to a two-dimensional elec- 
tronic subspace (it being assumed that other electronic states of the molecule are too far 
away to have any effect) and small displacements of the nuclear coordinates from some stan- 
dard configuration (so that only linear terms in the nuclear displacement coordinates appear 



in the Hamiltonian below). The solution to the mathematical problem (the DJTE) embod- 
ies the correlated nuclear-electronic trajectory near a conical intersection of the (diabatic) 
potential surfaces. Under these circumstances the usual Born-Oppenheimer approximation 
breaks down and the description of the combined dynamics is non-trivial. 

The total Hamiltonian consists of Hmoi for the internal degrees of freedom of the molecule 
and an interaction term with the environment Henv 

H = Hmol + Henv (13) 

The first term is a function of the electronic and nuclear coordinates, while the second 
term may contain also an externally imposed time dependence. Our restriction to a two- 
dimensional electronic subspace |1 >, |2 > removes from the formalism the presence of 
electronic coordinates and leaves only the nuclear coordinates. Two of these, designated 
ga,gb, are of interest. Hmoi when expressed in terms of the bosonic creation {a^ ,b^) and 
annihilation (a, h) operators of the nuclear motion, takes the form 

Hmol = |{a+a + b+b - ^[(«^ + a)^. " (b^ + b)cr,]} (14) 

Here u is the frequency of oscillation of the nuclear motion, and k is the electron-nuclear 
coupling strength expressed in dimensionless units. The 2x2 matrices ax and az are the 
familiar Pauli operators acting on the electronic |1 >, |1 > subspace. Equivalent represen- 
tations of the Hamiltonian Hmoi are given in works on the Jahn- Teller effect ( J4|, Isl, Isf ) ; i.e., 
in terms of the nuclear coordinates g^, q^ or of the associated cylindrical coordinates (g, 0), 
where Qa = qcoscj), qb = qsincf), as, e.g., in Eq. 3.5 of y. 



The algebraic expression for the ground state doublet proposed in J9|, |10|, and which 
solves the time-independent Schrodinger equation for the Hamiltonian Hmoi to a good ap- 
proximation, has the following (unnormalized) form: 

i^iQa, qb) = exp -- [(g„ - ka^f + {qb + ka^f] (15) 

(Intuitively, this form is suggested by analogy with the ground state solutions of displaced 
harmonic oscillators, but its justification is in its close agreement with exactly computed 
eigenvalues of Hmoi ^2.])- To obtain the ground state doublet we operate with ipiqayqb) on 
any two linearly independent combinations of the basic vectors |1 >, |2 >. In a column 



vector representation these are just ( j , (',')• The exponential, which includes non- 
commuting matrices, can be manipulated by use of the commutation relations between the 
Pauli-matrices to give, in terms of the cylindrical coordinates defined above, the expression 

'^^{q.t'P) = exp[— /c^ — q'^ /2][cosh.{kq)I — sinh(A;g) (ctz cos — ax sin 0)] (16) 



I is the unit 2x2 matrix. One notes that this is a single valued function of (there are no 
cos 0/2 terms), as indeed is required of the wave function of a closed system [3]. 

■ J ^\ 

By operating with ^ on ^ one gets the two degenerate ground- state functions 

\E'(g,0)(= \E'_) and \E'*(g, 0)(= \1^+). The eigenvalues and other related properties of these 
states have been calculated in |j|, [ij . Here we compute the phase-change for each function, 
as the angular coordinate changes by a full period between = and = 27r. 

One procedure to induce such change in an internal coordinate (and physically, perhaps, 
the only consistent one) is to consider it being guid ed by an external force along a circle. (The 



concept of a guiding potential was used in, e.g, jlj], but here we guide the angular coordinate. 



rather than the radial one.) To achieve this, one needs an external, time dependent agent 
acting on the otherwise closed system and this is the role played by H^nv in the Hamiltonian 
shown in equation (jlHj) . We suppose that there is an environment Hamiltonian which 
induces a de/ta-function like behavior in the wave function (forcing to equal Qt) and that 
this time-dependent H^nv dominates the kinetic energy of the angular variable. In this, 
delta-function limit the variable turns into a (classical) parameter and is no longer a 
" degree of freedom" . 

We thus get from equation (|T^ for the expectation value of the phase change in the \E'~ 
state the expression 

<[S]l'> = [ImJ dt[j dqq^t{q,nt)-^4q,nt) 

POO 

+ 5K]/ dqq\<il4q,nt)\'] (17) 

Jo 

where 6K is the second term on the right hand side of equation (J12j) involving the space 

differentials. In all integrals the integration is over the radial coordinate q only, since (p = Qt 

is treated as a parameter. 

It can be shown that 6K is identically zero for the both \l/_ and \E'+. (Remember that the 

gradient operator in 6K involves now only the radial degree of freedom q.) The evaluation 
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FIG. 1: Expectation value of the phase change after a full cycle in the coordinate space versus the 
coupling strength. We plot [< S >]q^ ( equation (fT7|) , with 5K = ) against k (introduced in 
equation ((T^ ). For large k, the phase approaches vr. 

of the first integral in equation (fTTj) leads to the plot shown in the following figure for the 
mean phase after a full cyclic revolution < [S]q^ > as function of the coupling strength k. 
(Fig. d) 

As seen in the figure, the acquired mean phase for \E'_ increases monotonically with the 
coupling and levels off for strong coupling {k » 1) to vr (the value of the Berry-phase). The 
corresponding phase for the partner state \E'+ is the negative of this value, and any linear 
combination of the ground state doublet will result in intermediate values between the two 
extremes ±7r. The phase depends only on the strength of the coupling k. It is independent of 
the adiabatic parameter ^, since the integrand contains only the instantaneous value of the 
initial component \I'_ , and no admixture from its partner \E'+ . (Applying our equation ()17|l 
to the eigenstates ^^ of yy] expressed in a coordinate representation reproduces exactly the 
results obtained in that paper. However, evaluation of the expectation value of the phase- 
shift operator proposed in pj for the states \E'^(g,0) in this work, where the rotating- wave 
approximation is not made, yields values that diverge quadratically for large k.) 

We conclude with an interpretation of the " closed-system" phase. In this we follow p, |S| . 
For low values of the coupling constant, the wave function is smeared over the origin g = 
and cannot be said to circle around this point, which is a point of degeneracy of the two 
states. Then there is hardly any acquired phase. For large values of the coupling, the wave 



function is located near q = k, meaning that it keeps away from origin, so that circhng 
around it can achieve the full measure of the geometric phase. 

On the other hand, it has been known for some time that (in the adiabatic limit) the 
phase change comes about abruptly, precisely at the moment of circling when a component 
amplitude vanishes. (This occurs when cos^ = 0, or Qt = n. The abrupt change is clearly 



seen in the figures of |l6l.ll7| and has recently formed the subject of a Letter |18|.) By the 
interpretation just given, the phase change is a measure of the extent that a circling in the 
coordinate space scans the zeros of the wave function in the region encircled. Since zeros 
(nodes) in the wave function are known to affect (in general, raise) the energies of the states, 
it is natural to find that the phase acquired during a revolution determines the ordering of 
the energy levels. Such connections between phase change and energy levels have been noted 
first in fTJ and more recently in Iq]. 
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